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INVESTIGATIONS WITH SATELLITE DATA 
TEMPERATUBE RETRIEVALS 


Introduction 

It can be shown that there is no unique solution (e . g. Fritz, 1969) 
to the problem of atmospheric temperature retrievals from satellite radiance 
measurements . The problem which does remain is to find the "best" solution 
from the infinite set that exists. 

For this purpose there is some merit in beginning with simulated "data" 
instead of real data, because one of the major difficulties with real data 
is that the true answer is not precisely known; real data contain various 
errors and uncertainties themselves. By contrast, for simulated "data", 
the true answer is known. And, further, all physical quantities which 
enter into the problem can also be simulated without error, a fact which 
is never true about the real atmosphere. 

Because of these attributes, the processes used to retrieve the atmos- 
pheric vertical temperature distribution can be studied with simulated "data". 
And that process which produces temperatures closest to the known true tempera 
tures can be judged to be the best from among those processes tried. 

Almost all retrieval methods involve the multiplication of the measured 
radiances by some numbers. In particular, linear methods usually involve 
the search for a set of coefficients, F(k,j) , such that 

AB(k) » Z F(k,j) -AR(j) 

j 

will achieve the "best" solution. See Appendix (A) . B(k) = Planck function 
at the pressure level, k, at a reference frequency, usually taken at 700 cm 
when the 15 ym CO 2 band is used in the satellite measurements, 

R(j) is a radiance at the frequency, j. 
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The subscripts denote the following: 

"g" = a guess (or initial estimate) 
m = the measured value 
r =* the retrieved value. 

Because B(k) depends only on the temperature, T(k) , it is easy to 
compute T(k) from the inverse of the Planck function, once B(k) is known. 

Two linear methods, the "minimum information" method and the regression 
method have been used operationally. The regression method defines the "best" 
coefficients, F(k,j) as the set which produces the smallest temperature error, 
in a least square sense, in a set of radiosonde data selected from a wide 
area during some time interval before the current satellite data are obtained 
at a particular place. A similar objective, but with additional assumptions, 
to make the problem more tractable, is employed in the minimum information 
method. In view of these limitations should one expect that the set of 
coefficients will produce the "best" temperature retrieval at a particular 
time and place? A recent paper (Fritz, 1976, see Appendix A) suggests that 
better results can be obtained. This is achieved by using the most recent 
radiosondes available and modifying the coefficients F(k,j) to a new set, 
G(k,j) so that retrieved temperatures agree exactly with one radiosonde 
temperature profile. Thus, instead of trying to "beat" or even match the 
radiosondes. We "join" the radiosondes and use the satellite radiances to 
improve temperature determination between the radiosonde stations and also 
in time until the next radiosonde launch times. This final report expands 
on a previous report which is included here as Appendix A. Appendix A, 
which has been submitted to the Journal of Applied Meteorology (JAM) , sum- 
marizes the results obtained from simulated "soundings" for "stations" 
between Omaha, Nebraska and Springfield, Illinois. In this final report, 
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additional results are discussed; it also includes some findings obtained 
after the earlier report was submitted to the JAM. 

Data Used 

To investigate the effect of the coefficients two sets of "data" were 
used. One set depended on ship temperature soundings in the GATE project. 
Another depended on temperatures between Omaha, Nebraska, and Springfield, 
Illinois. The Omaha profile (which served as the "guess" profile in the 
retrieval method) and the Springfield temperature profile are sho\m in 
Fig. 1. 

Effect of Noise Parameter on Minimum-Information Retrievals 

With the Omaha temperatures as an aid , and with the transmit tan ces » and 
their corrections for temperature supplied by Dr. M. Weinreb (NESS/NOAA) , 
the coefficients F(k,j) used in the minimum information method, were 
computed. The F(k,j) depend on the noise parameter, included in the 

retrieval (See Appendix A eq. 4). 

Fig. 2 is a plot of the average temperature error, AT averaged over 
the various layers indicated, as a function of the for the station pair, 
Omaha (first guess) and Springfield, Illinois (retrieval). In deriving 
Fig. 2, the "measured" radiances contained no errors. We note that as 

approached zero, AT often became large. Thus for the layer 1000-800 mbs, 
AT = 10. 4K; for the layer 800-600 mb, AT was -13. 3K. And, for a large, AT 
approaches the difference between the first guess temperature and the true 
temperature, because F(k,j) and AB(k) approach zero (see eq. 4 in Appen- 
dix A). At least for the one sounding (Springfield) if we accept the 
interpolations in Fig. 2, there was no single value of cr^ which would 
make AT = 0 for all layers. For the layer 1000-800 mbs apparently 
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a =1.0 would have made A T = 0, for the layer 600-400 mbs, a - 0.15 

N 

would have been best, while for the other layers there may not have been 
any value of which would have made AT = 0. Thus since F(k,j) is 
sensitive to the value of a , it seems evident that it would be impossible 

iN 

to find a set of F(k,j) which could '.make AT = 0 for all layers. In view 
of this defect, it seems advisable to seek a set of coefficients, G(k,j) 
which would make AT = 0, at least at Springfield. Apparently such a set 
of G(k,j) cannot be obtained by adjusting and possibly not even by 
adjusting other parameters in the minimum information method. 

Vertical Distribution of Coefficients 

Since we employed the six 15 pm CO 2 channels used in the VTPR instru- 
ments, at each pressure level, k, in the atmosphere, there are six coeffi- 
cients, one from each spectral interval, j; i.e., at each pressure level 
there are six values of F(k.j). 

There are an infinite number of ways in which the F(k,j) can be 
modified to make AT = 0 at Springfield. A few were tried so far in this 
study. 

These included the following methods: (see Appendix B, for the equations 
for the various methods . ) 

(1) Change one of the F(k,j) so that AT = 0. To do this the largest 
(in absolute value) of the six F(k,j) at each level was changed. 

(2) Change all the F(k,j) by-multiplying them by a constant; thus 

G(k,j) ^ C(k).F(k,j) 

This method used in Appendix A. 

Methods 3 and 4. 

Make the Sum of G(k,j) = 1.0. Tnis was done in two ways. 



(3) Change the two values of F(k,j) for which the absolute values 
of F(k,j) were the largest. As in Methods (1) and (2) above, 
certain instabilities occur in the resulting G(k,j). There- 
fore Method (4) was devised. 

(4) Change the highest absolute value of F(k,j) and also the one 

for which AR(j) was the furthest away from ’ 

lARCj^ax^ -4R(j)l = maximum 

This makes the G(k,j) stable with k. 

The vertical distribution of F(k,j) (minimum-information) and G(k,j) 
(adjusted) are shown for several spectral channels and for several methods 
in Fig. 3 

The coefficients, F(k,j) for three spectral channels are plotted in 
Fig. 3(a) against the vertical coordinate, k. (The other channels are not 
shown to avoid complicating the diagram) . These are the coefficients ob- 
tained in the minimum-information retrieval method, when the temperatures 
at Omaha are used as a first guess in order to apply temperature corrections 
to the transmit tan ces ; the vertical variation of the transmittances , to a 
large extent, determine the coefficients. 

For channel 6, v =747.6 cm~^, the largest coefficients are found 
near the ground. Thus at p =.1000 mbs, F(100,6) = 4.9. We shall find 
below that this is a large value, and is doubtless responsible for the 
large temperature errors produced in the minimum- in format ion method. The 
other channels shown have little influence near the ground: however, chan- 
nel 5 (not shown) has an influence almost as great as channel 6. 
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Near the 150 mb level, channel 3, v = 695.3 cm ^ , has the greatest 
influence of those channels shown; and finally, channel 1, v = 668,5 cm 
has the greatest influence near the top of the atmosphere , near the 2 mb 
level . 

By comparison, the coefficients G(k,j) are shown in Figs. 3(b), 3(c) 
and 3(d) for several other retrieval methods. 

Fig. 3(b) shows the results for channel 6. In Fig. 3(b) curve 3 is 
the same as F(k,6) in Fig. 3(a). This is the minimum -in formation method 
curve. The coefficients G(k,j) derived from 3 of the adjustment methods 
are also shown in Fig. 3(b). The significant result is that in the lower 
atmosphere, near p = 1000, the adjustment coefficients are much smaller 
than F(k,j). The values of G(100,6) vary from about 3 down to about 0.5, 
whereas F(100,6) is 4.9. This indicates that in order to obtain small 
temperature errors below the 900 level, the minimum- information coeffi- 
cients must be substantially reduced in the lower atmospheres. 

Fig. 3(b) also shows the instability in the G(k,j) introduced when 
the denominator in the adjustment computations get too small (see Appen- 
dix B) . Thus for curve number 4, near the 270 mb level, the coefficients, 
oscillate rapidly between about -3 and +8.1. However this oscillation is 
eliminated in Method 4 as shown in curve B, for which EG (k,j) = 1. Fig. 3(c) 
shows similar results for channel 3, v = 695.3 cm . Here curve (C) is 
the same as F(k,j) in Fig. 3(a) for channel 3. Fig. 3(c) shows how un- 
stable some of the adjustment methods can become. By contrast, curve B, 
which shoxjs the G(k,j) for Method 4, shows how stable the coefficients 
are. Still the relatively small differences from the minimum-information 
coefficients are required to reduce the temperature error to zero at 


Springfield. 
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Finally Fig. 3(d) shows similar curves for channel 1, v = 668.5 cm 
Here, of course, the major changes occur near the top of the atmosphere. 
Curve A, adjustment method 3, gave poorer results; the coefficients are 
large, because the denominator in the calculation of G(k,j) are too small 
in Method 3 for chat part of the atmosphere. Again Method 4, yielded 
stable coefficients and gave relatively small errors (see Figs. 4 and 5). 
Making EG(k,j) = 1.0 can be justified some^^hat by the following 

j 

reasoning: 

Consider an isothermal atmosphere; then 

R (j) = SB (k,j)AT(k,j) = B (k,j)SAT(k,j) = B (k,j) = B (T ,j) 
t fct tic ^ 

since SAx(k,j) = 1, and T is independent of height (subscript 
k t 

"t” = true). 

If we select also an isothermal first guess profile, B (k) , then, 

O 

similarly, R (j) = B (T ,j), and 
S o o 

AR(j) = R.(j) - R^(j) = B (T ,j) - B (T ,j) = B (T ,700) - B (T ,700) 
t g ct S8. tc gg 

= constant 

Therefore AR(j) is nearly constant and nearly independent of frequency, 

j . 

Moreover, B(k, 700) - B_(k,700) = AB (k,700) = ZG(k,j) *AR(j) 

C a ^ 

Therefore AR(j) = AR(j) SG(k,j) 
and SG(k,j) = 1. 
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Similarly, ZG(k,j) would doubtless also equal about 1.0 if the first- 
guess profile differed at all levels by a constant temperature difference, 
say, 3K, from the true temperature for a realistic non -isothermal atmos- 
phere. The reasoning would be similar to that for an isothermal atmos- 
phere; but the answer would not be as accurate because the Planck Function, 
B(k,j) is not exactly linear with temperature. If the first-guess tempera- 
ture profile differs in arbitrary ways from the true profile, the assumption 
of SG(k,j) = 1 might be less justified. 

It should be emphasized that for all these methods, AT = 0 , exactly, 
at Springfield. However the application of the coefficients, at "stations" 
between Omaha and Springfield, did not yield identical results because 
the AR(j) were, of course, different for the intermediate stations from 
their values for Springfield. The results are given in Table 1. 

Table 1 shows that the average absolute error of the "adjustment" 
retrieval temperatures averaged over 200 mb layers is smaller than the 
minimum-information retrievals at all levels except the 200-5 mb strato- 
spheric layer, x^here for most methods the error was small anyway. Most 
striking, however, is the improvement of all four adjustment methods, over 
the minimum-information method, in the layer from 1000-800 mbs. Whereas, 
the minimum— information method makes essentially no improvement over the 
first-guess error of 3K, the adjustment method Nos. (2) and (3) reduce 
the error to less than 0.5K. Method (2) involves division by 

[B (k) - B (k) 1 = AB . This quantity approaches zero xfhen the retrieved 
r g rg 

temperature is approximately equal to the first guess temperature. Under 

those conditions, AT^ may be large and oscillate in sign. Yet as can be 

seen from Table 1, AT averaged over 200 mb layers was small. This is a 

result of the fact that when B (k) = B (k) the sign of (B ™ B ) Xus elf 

^ S v g' 
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TABLE 1. Average Absolute Temp Error (°K) 



Adjustment Method 

Ap 

(mb) 

FG- 

True 

Min Info 

■ 

(Best) 

2 

3 

4 

.000-800 

3.0 

2.9 

[1.0] 

0.3 



800-60U 

7.4 

2.3 

[1.6] 

1.2 



600-400 

MEM 

1.6 

[1.9] 

1.0 


1.0 

400-200 


3.8 

[1.8] 

1.2 


1.3 

200-5 

1^1 

0.1 

[0.1] 

0.2 

0.8 

0.2 


Adjustment Method (see p. 4, 5) 
*[(1), G(k,j) = F(k,j); G(k,j = C*(k).F(k,j ] 

(2) , G(k,j) = C(k)*F(k,j) for all j 

(3) , SG(k,j) = 1.0; G(k,j) = F(k,j); G(k,j ; G(k,j 

(4) , 2G(k,j) = 1.0; G(k,j) = F(k,j), ni^); G(k,j) 

lAR(j J) - AR(j’)j = maximum 
*[for Method No. 1 

(Omaha + 90 km) was used for the adjustment instead of 
Springfield, Illinois. A test showed that even when this 
station was used with all F(k,j) ch^ged, the result was more 
similar to the results shown above for Method (2). Also, it 
should be noted that in that early test, Omaha, the first 
guess station, was inadvertently made super adiabatic from 
the surface to 700 mb; the temperature at 700 mbs was about 
9K colder than in the later tests.] 
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oscillates in sign witii height about the critical height. Thus in Fig. 1 

at about 275 mbs, B^(k) = B„(k) ; but just above that level AB is negative 

r g 

and just below 275 mbs AB^^ is positive. Therefore, if one averages over 
the layer from 400-200 mbs, AT appears to be quite stable. However for 
individual pressure levels, which should not in any case be used for 
temperature retrievals, AT can be large. For example. 


Temperature Error (°K) 


Station 

P(mb) 

Omaha 
+90 km 

+180 

+270 

+360 

+450 

285 

+3.5 . 

+2.4 

+2.7 

+3.7 

1.0 

271 

+25.1 

+1.9 

+4.6 

+22.0 

+5 .2 

258 

-11.0 

+2.6 

+2.9 

-8.4 

-1.8 

245 

-5.4 

+2.5 

+3.7 

-3.9 

-1.0 


At 271 mbs, AT = 25.1 at station Omaha + 90 km, and 22.0 at station 
Omaha + 360 km. Clearly near the 271 mb level, the retrievals are sensi- 
tive to the radiances because the G(k,j) have large absolute values. 

To overcome such large AT oscillation at particular pressure levels, 
we may average the values of G(k,j) over a few levels. But it seemed to 
be preferable to use Ifethod (4) mentioned under Table 1. With that method 
we get more stable values, as for example 


Temperature Error (°K) 


Station 

mi 

Omaha 
+90 km 

+180 

+270 

+360 

+450 

285 

-1.1 

2.8 

2.5 


0.7 

271 

-1.0 

2.8 

3.0 


0.7 

258 

-0.9 

2.7 

3.4 

0.1 

0.6 

245 

-1.0 

2.7 

3.8 ' 

0.2 

0.7 
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The GATE Experiment 

Before the Omaha-Springfield experiment was run, an experiment was 
run involving soundings from the Canadian Ship QUADRA, located at about 
9.3N, 22W, one of the ships in the GARP Atlantic Tropical Experiment 
(GATE). Atmospheric temperature soundings were available for'lSOOZ on 
eight days beween the period July 5, 1974 and August 8, 1974. For the 
eight days, data for one day was used as the first guess, and data for the 
next available date were used as the truth [actually the order was inadvert 
ently reversed so that, e. g., data for July 8 were used as the first guess 
for retrieval of July 5 data; the computation was not changed, because for 
the purpose of this study, the order does not make any difference.] In 
this early experiment, comparison was again made between the retrievals 
for the minimum-inf ormation method, and the adjustment method. For the 
adjustment method the coefficients F(lc,j) were modified by two of the 
different ways discussed above; namely. Methods (1) and (2). The results 
are sho\m in Table 2. 


TABLE 2 

Average Absolute Temperature Error (°K) 
for GATE Ship Quadra 



Method 

Adjustment 

AP(mbs) 

FG-True 

Min-Info 

1 

2 

lGOO-300 

1.8 

3.3 

0.5 

0.5 

800-600 

1.3 

2.2 

0.8 

0,7 

600-400 

0.4 

0,4 

0.6 

0.7 

400-200 

0.6 

0.9 

1.2 

0.8 

200-5 

2.8 

0.6 

0.4 

0.4 

FG = First Guess 

Min-Info = Minimum-information method 
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In Table 2, again, the minimum-information method yields poorer re- 
sults than even the first guess in the lower layers, below 600 mbs. But 
the adjustment method improves over the first guess. (It should be noted, 
however, that the ship’s observations might have been in error; on some 
days the average temperature in the layer 1000-800 mbs, differed from the 
adjacent day’s temperature by over 3.5K; this seems like a rather large 
change in the tropics. Since the radiances in this study were simulated, 
any measurement errors on the ship would not affect the results. In a 
real case, however, ship temperature measurement errors would lead to 
inconsistencies with satellite observed radiances, and apparent "errors" 
would show up in retrieval comparisons with radiosonde temperature measure- 
ments) . 

intermediate levels, between 600 mbs and 200 mbs, the fxrst guess 
error (i.e., the observed variation from day to day) is so small, namely 
about 0.5K, that none of the methods can improve on the first guess. In 
this pressure interval, the minimum— information method did as well, or even 
better than the adjustment method. 

For the layer from 200-5 mbs, the first guess error was rather large. 
This was caused by a rather poor add-on of temperature profile for the high 
layers above the radiosonde level* In this layer the retrievals by all 
methods were rather good. This is a bit surprising, perhaps, because this 
layer contained the tropopause, which was located at about the 100 mb level. 
Usually, the minimum- information method yields large errors near the tropo— 
pause. Perhaps, in this study, the layer 200-5 mbs is large enough to 
eliminate errors which appear in thin layers, or at particular pressure 
levels . 
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Results for Other Layer Thicknesses (Omaha-Springf ield) 

There is considerable discussion in the literature regarding the 
vertical resolution, or vertical resolving power, of the temperature 
retrievals from satellite data [e.g., Conrath, 1972; Thompson, et , 1976], 
Everyone concludes that there is a trade-off between the vertical thickness 
over which the retrievals are averaged, and the error of the temperature 
retrievals. In this study the temperature averaged in the vertical over 
various pressure layer intervals was computed. The results are shown in 
Fig. 4. The figure shows the average absolute temperature error averaged 
over the ftve stations * between Omaha and Springfield, for various layers 
above lOQO mbs , i. e., the ordinate 900 mbs means the layer from 1000-900 mbs, 
the ordinate 800 mbs means the layer from 1000—800 mbs, etc. Figure 4 shows 
that for layers thinner than 1000-500 mbs, the minimum-information results 
were substantially poorer than for any of the adjustment methods. Actually 
for the layer from 1000-900 mbs, the minimum-information retrievals were 
much poorer than the first guess error; for, whereas the first guess tem- 
perature error was about 1.5K, the minimum-information retrieval was more 
than 4K for that layer. Even for the layer 1000—800 mbs (as already noted) 
the minimum-information retrievals did not improve over the first guess. 

It is not until the thickness reaches 1000-500 mbs that the minimum-infor- 
mation retrieval errors are approximately as small as the adjustment method 
retxievals. This suggests that in the lower atmosphere, the vertical reso- 
lution, with the minimum-information method, for average absolute tempera- 
ture errors of about IK, is about 500 mbs. 

However, Fig. 4 shows that with adjusted coefficients, the temperature 
error is less than IK even for layers as thin as the 1000-900 mb layer. 
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From the theoretical results which estimate the vertical resolution on the 
basis of the transmittances and the 'Veighting functions (Conrath, 1972) , 
it is not likely that the good results for the adjustment method depends 
on any inherent resolving pox^er in the satellite radiances. Rather, the 
adjustment method is doubtless taking advantage of the similarity in the 
structure of the vertical temperature profiles at the interpolated sta~ 
tions" and at the nearby radiosonde stations. 

For the average absolute temperature error for the layers thicker 
than 1000-300 mbs, the results seem more erratic. We note that Method 3 
(ZG = 1) shows the poorest results for the layer 1000-100 mbs. This 
happened because the coefficients G(k,j) became very large in the layers 
involving channels 2 and 3. For those channels, the AR were respectively 
-10.33 and -10.56; the difference between these is the small value 0.228. 
Since this small number, 0.23, appears in the denominator at those levels, 
k, x^hera F(k,j) are maximum for channel 2 and next to the maximxmi for 
channel 3, the coefficient G(k,j) will be large in the height intervals 
where channel 2 and channel 3 dominate, i^^hen ; those large coefficients, 
derived from the Omaha-Springfield pair, are then applied to other stations, 
the results are likely to shoxv large errors. Moreover in this method of 
finding G(k,j) unlike Method 2, for which G(k,j) = C(k)*F(k,j), the large 
values of G(k,j) will persist over a deep layer as long as F(k,j) is the 
largest of the six values at any level, k. Thus, for example, in the 
"Omaha + 90 km" sounding, the temperature error x^as between -4.9K and 
-11.7514 at all levels between 233 mbs and 141 mbs. Such large tempera- 
tures of the same sign for such a deep layer xd.ll produce large temperature 
errors even if averaged over a deep layer, and will also adversely affect 
the height errors as we shall see presently. By contrast, large temperature 
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errors, of even 20K, in Method 2, [G(k,j) = C(k) F(K,j)], do not produce 
poor results over layers, because the sign of the error for adjacent 
values of "k'V generally oscillate in sign. 


The Height Error for Various Pressure Levels 
In this study, the height was assumed known at the 850 mb level, 
where the height was taken to be 1330 m. The height Z(k) for various 
pressure surfaces, k, was computed from the 850 mb height and the re- 
trrlevad temperature structure from the formula 

Z(k) = Z(k-l) - (| ) [T(k) ■¥ TCk-l) . ! 

* 

— = 29,287 
g 


Tlie results, in Fig. 5, show that Method 4, [SG(k,j) = 1, MOD], and 
Method 2 [G(k,j) = C(k) *F(k, j) ] gave the best results. At all pressure 
levels up to 10 mbs, the average absolute height error did not exceed 20 m; 
in general, the second of these two methods gave somewhat better results 
than the first method. The other methods illustrated in Fig. 5, show that 
the minimum-information method was poorest at 500 mbs although all methods 
improved substantially over the first-guess height difference from the true 
heights. At the 100 mb level, Method 3, [2G(k,j) = 1] (t^'^o largest coeffi- 
cients changed) gave the poorest, results, and even did not improve over 
the first guess height error. This resulted from the fact that, as stated 
earlier, large temperature errors, of the same sign persisted over relative 
ly thick layers. 

From the results in Figs. 4 and 5 and from Table 1, it appears that 
i-Iethods (2) and (4) gave the best results; and if, for some reason, it is 
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necessary to avoid large oscillations in the temperature for individual 
pressure levels. Method (4) would be the one recommended. 

Values of F(k,.j) AR(j) and B(k.700) - Their Relationship 

Linear Methods 

As stated earlier, linear methods generally solve the satellite tempera- 
ture retrieval problem by solving for 

B (k) = B (k) + AB^(k) 

^ s 

AB = SF(k,j) -AR(j) 

^ j 

AR(j) = \(j) - Rg(j) 

In the minimum-information method, an attempt is often made to make 

B (k) as close as possible to the true value. This can be done, for ex- 
§ 

ample, by selecting NMC’s latest 12 hour forecast at the place where the 

satellite measurements are being made. This would also have the tendency 

to make AR(j) small, because R^(j) would approximate Rj^(j)* particular 

if all B^(k) = Bj.(k) (although this would not be known to the experimenter), 

AR(j)- 0 and it X'Tould not make any difference what values F(k,j) had. [Even 

if B (k) = B (k) , AR would not equal zero exactly, because, the T(k,j) used 
g t 

to compute R (i) would not be exactly the same as the value in the real atmos 

g 

phere; the real t ' s are the only ones which the satellite experiences when it 
measures R^(j) ]. In practice, of course B^Ck) does not equal B^(k) , and 
AR(j) 3 ^ 0. Therefore, the values selected for F(k,j) become important; the 
set of F(k,j) selected will strongly influence the resulting values of B^(k) 
and therefore of T^(k) . Since a value of AR(j) = 0 would most likely result 
from a value of Bg(k) - B^(k), the retrieval error would tend to be smaller, 
AR(j) is . 


the smaller 
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However in the adjustment method in ^diich, for example, 

G(lc,j) * C(k)-F(k,j), the retrieval error does not necessarily get 
smaller as AR gets smaller. Take, for example the results for the 
Omaha-Springfield experiment reported above. The retrieval error is 
e.cactly zero in two instances; once when AR is a maximum, i. e., at 
Springfield, and again when AR = 0, and also (B - B ) = 0, i. e., at 
Omaha. The error is zero at Springfield only because the G(k,j) have 
been so selected;, so that for the given values of B (k) and therefore of 

9 

O 

snd of AR(j ) , (all of which are not optimum from the minimum— inf or~ 

mation method point of view), B^(k) = B^. (k) . ; Now as we proceed from- 

• ^ 

Springfield towards Omaha, AR(j) generally becomes smaller and B (k) 

s 

approaches Bj.(k) for the particular ’’station". Still, the adjustment 
method derives better retrievals, probably because the influence of the 
G(k,j) is dominant over the AR(j) and B^(k). However, when we arrive at 

O 

the station Omaha + 90 km" the minimum— information retrieval is better 
than the adjustment” method retrieval for deep layers, as indicated by 
the 510 mb height error (in Table 3 of the Appendix A) . At the "Omaha 
+ 90 km station the height error for the minimum— information method is 
14 m, for the adjustment method it is 23 m. At all other stations further 
from Omaha, the height error is lairger for the minimum-information method. 

In summary there are two main factors which affect the retrieval 
accuracy with simulated data; namely B and AB . And AB in turn depends 

*Table (3) shows the values of AR(j) and of N(j) for each station between 
Omaha and Springfield. Although the|AR(j)l generally decrease from Spring- 
field to Omaha, they do not do so nionOtomically ; this is doubtless due in 
to the random nature of the noise added to the radiances, and partly 
to the fact that the gradient of temperature was toward Omaha (low temp)* 
in midtroposphere, but was reversed above the 200 mb level. The radiances 
are affected by temperatures over a deep layer. 



Chan. 

No . 

Springfield 

Omaha 

1 

- 6.53 

2 

- 10.33 

3 

- 10.56 

4 

- 1.63 ' 

5 

2.37 

6 

3.29 


Springfield 

.11 

-.04 

.18 

.21 

.05 

-.17 


FG Spring- 

Omaha field 

71.5 65.2 

66.8 56.4 

63.1 52.5 

66.6 65.0 

89.5 91.9 

110.2 113.5 


AR(j) = CR - R^) 
; m g 


Station 


Omaha + 
90 km 

- 1.30 

- 1.69 

- 2.18 

- 0.46 

- 0.79 

2.36 


+180 

+270 

+360 

• 1.98 

- 2.57 

- 6.23 

- 3.90 

- 4.21 

- 8.64 

- 4.22 

- 4.36 

- 9.18 

- 0.10 

- 0.04 

- 1.74 


2.97 1.88 

3.16 3.23 


- 5.67 

- 8.27 

- 8.86 

- 0.90 

2.61 

3.68 


Radiance 

Noise 




-.09 

.43 

-.15 

-.73 

-.16 

.22 

-.08 

-.37 

-.08 

.31 

-.08 

.16 

.12 

-.06 

.26 

-.10 

.41 

-.06 

-.17 

.49 

-.28 

-.04 

.02 

-.47 

.13 

.21 

.47 

-.04 

-.10 

.28 

Radiances 




70.4 

69.8 

69.2 

65.5 

66.1 

65.1 

62.9 

62.5 

58.1 

58.5 

60.9 

58.9 

58.7 

53.9 

54.2 

66.2 

66.7 

66.7 

64.9 

65.7 

90.3 

91.8 

92.5 

91.4 

92.1 

112.5 

113.4 

113.3 

113.4 

113.8 
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on F(k,j) and AR(j). Moreover AR(j) depends on R^ and therefore also 

Srdc i ) 

on B . The F(k,j) depend on — r — — , where x is an atmospheric pressure 
S o X 

parameter. And since the t’s are temperature dependent, and in practice 
are adjusted with the first guess temperatures, even the F(k,j) also depend 

somewhat on B . Therefore B enters the problem in three wavs in the mini- 

g g • • 

mum-information method. And since in the Omaha-Springfield set, the Omaha 
first guess is furthest from the true value at Springfield, the errors in 
the minimum-information method might be largest at Springfield when Omaha 


is used as the first guess, 


However, in the adjustment method used in this report, the temperature 
error is forced to be zero at Springfield, just where the first-guess error 
is largest. And since the vertical temperature profile structure is pro- 
bably more nearly like the Springfield profile for the "stations'* nearest 
to Springfield, one might expect that the coefficients G(k,j) would have 
an opposing effect to the effect of B . Thus the computed temperatures 

O 

retrieved from the adjustment method, at a particular "station" between 
Omaha and Springfield, would depend on the relative influence of (B^ - B^.) 
at the station vs. the influence of the G(k,j) which have their maximum 
effect at Springfield (where Omaha is used as the first guess). In the 
test performed in this report, judging by the 510 mb height error com- 
parisons in Table 3 of Appendix A, the influence of the G(k,j) seems to 
have dominated except at the station "Omaha + 90 km". Perhaps at that 
station the smallness of (B^ - B^) might have dominated. However, more 
tests will be needed before the i argijments presented here can be accepted. 
Preliminary results, based on a reversal of the roles of Springfield and 
Omaha, raise some questions about the relative roles of B (k) and G(k,j). 

g 
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An indication of the composite influence of the various factors is 
shown in Fig. 6. This figure shows the temperature error for each station 
for layers 200 mbs thick, namely, 1000-800 mbs, 800-600 mbs, 600-400 mbs, 
and 400-200 mbs. The station numbers refer to the distance from Omaha, 

e. g., number 2 is 90 km from Omaha, number 3 is 180 km from Omaha, etc. 
Figure 6 shows that the temperature errors near stations 5 and 6 were 
always small for the adjustment method. This is undoubtedly the in- 
fluence of the G(k,j) which had been selected, so that with Omaha as 
the first guess, the retrieval error would be zero at Springfield (90 km 
further from Omaha than station 6). The minimum- information method, 
almost always, gives a larger temperature error than the adjustment 
method especially at stations 5 and 6. As we approach Omaha, however 
the minimum- in format ion method error becomes smaller, relative to the 
adjustment method. In fact, at station No. 2, in two layers, namely 
the layers 800-600 mbs and 400-200 mbs the minimum-information method 
error is somewhat smaller than the adjustment error. At station No. 2, 
the first guess error is smallest, AR(j) is smallest, and the influence 
of the F(k,j) and of the G(k,j) is reduced since they are multiplied by 
the AR(j) to yield AB'^(k) . 

Interestingly, the retrievals in both methods, yield either maxima 
or minima at station No. 4, midway bet^<?een Omaha and Springfield. The 
precise reason for that behavior is not obvious. Doubtless it repre- 
sents the interacting influence of the three factors, B^(k) , AR(j) 
and F(k,j) or G(k,j) as the distance from Omaha increases. 
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Some Philosophical Questions 

In the linear methods of retrievals the method often involves finding 

Br(k) = BgCk) + AB(k) 


AB(k) = B^(k) - B^(k) = 2F(k,j) R(j) = ZF(k,j)[R (k) - R (k) ] 
r g -1 * *1 o 


The basic idea is to select, in some way, the coefficients F(k,j), which 

when multiplied by the AR(j) will give the "best" value of AB^(k) . This 

value of AB (k) is then to be added to B (k) so that B (k) will be as close 
r § ^ 

to the true value, B^(k) , as possible. 

However to a certain extent the coefficients F(k,j) are selected inde- 

T T “1 

pendently of AR(j) . In the minimum-information method F(k,j) = A [AA + yl] 
(see Appendix A) where 


A(k,j) 


Y 


= (k>.1 ) 

dx(k) 


= ( — ) (see eq . 4 in Appendix A) . 

^T 


"x" being a height parameter, usually some function of pressure, such as 

Since the t’s are temperature dependent, the A’s are also temperature 
dependent; the A's therefore depend on the first guess temperatures, which 
are generally not equal to the true temperature. Moreover in the real case 
the A(k,j), even for a standard atmosphere, may have been computed from 
theoretical evaluation of the CO 2 15 pm spectrum and may not be strictly 
correct; (but, in the simulated procedures, this should not be a factor). 

In addition, the selection of the value of y will strongly influence the 
value of the F(k,j). Normally, the "y" is chosen to be compatible with 
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the values of the noise in the satellite instrument. But we found in this 
research that even if the noise is zero in the radiances , and we set y = 0 
the temperatures near 1000 mb level in the minimum-information retrievals 
contain very large errors, because the F(k,j), with k corresponding to 
p = 1000 mb, have very large values* oscillating in sign with frequency, j. 

One may ask then, what value of y should be selected, when the noise 

2 

in the radiances is zero? In this research a value of y = (0.25) /lOO 
yielded what seemed to be stable values of AT^ of about 3 K near the 
surface. Values of y = 0 gave unrealistically large errors. And values 
of y very large make the retrievals approach the first guess. These re- 
sults were found when the radiances contained no error, and also when 
the radiances contained typical instrumental noise errors. Therefore, 
the value of y seems to be somewhat arbitrary and is not linked closely 
to the instrumental noise, except possibly when the instrumental noise 
is very large. This also, therefore, implies that the coefficients, 

F(k,j) are somewhat arbitrary and not necessarily optimum for a part- 
cular atmospheric situation. Therefore the method of adjustment proposed 
here, seems to have merit, from the philosophical point at least, since 
it adjusts the F(k,j) to radiosonde temperatures near the times and place 
of the satellite measurements. 
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Some values of F(100,j) for the Omaha sounding as function of 
are shown in the following table: 



1 1 

2 

3 

4 

5 

6 

0 

+.2444+01 

+.2399+01 

-.3757+01 

+.1968+02 

-.7772+02 

+.7734+02 

0.050, 5 X 
0.025 

+.1887+01 

+.1106+01 

-.1691+01 

+.1299+02 

-.5608+02 

1 +.5735+02 

0.50, 5 X 
0.25 

-.3996-01 

-.5162+00 

+.1333+01 

-.2492+01 

-.5643+00 

+.4794+01 

5.0, 5 X 
2.5 ; 

-.1548-02 

-.1087-01 

-.1930-01 

-.2968-01 

.1 

1 

+.1924+00 

i 

+.2624+00 

1 

Adjusted I 
Value j 

+.3092-02 

-.3994-01 

+.1031+00 

-.1928+00 




The first value under column gives the value for Channel 1. Then the 
five noise values used for the other five channels are given. The values 
of Ojj are the nvnnbers used in eq. 4 (Appendix A) and are not the values 
added to the computed true radiances. In this table the numbers +01, 

+02, -01, +00 refer to the exponent of 10; thus +.2444+01 means 2.444. 


We note that the coefficient in all channels increase in absolute magnitude in all 
channels with decreasing 0 ^^, when p - 1000 mbs. And when the large coefficients 
are applied to the radiances, large errors in temperature are computed. Since 
it would be difficult to select the "best" set of F(100,j) the adjustment 
method which yields G(100,j) seems like a reasonable approach. 

Hurricanes 

It would be nice if something could be done to measure temperature changes 
which occur around tropical disturbances which develop into hurricanes. If 
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the temperature changes were large enough and could be eventually monitored 
from geostationary satellites perhaps forecasters would have a tool for 
forecasting hurricane development. However, since the central part of 
hurricanes, except for the eye, is always occupied by extensive cloud 
cover, infra-red measurements cannot be expected to indicate much about 
the temperature changes inside the cloudy area. 

We therefore collected radiosonde temperatures for areas surrounding 
Hurricanes Camille (August 1969) and Celia (August 1970) in various stages 
of development. These data were plotted in several ways. For example at 
the 700 mb and 500 mb levels polar diagram charts were prepared which 
showed the temperature distribution at all azimuths and at distances from 
less than 200 km to about 1000 km from the storm center. In order to get 
sufficient coverage, the data were plotted for several days, on one chart, 
all relative to the storm centers. For Camille, this was done separately 
for the tropical disturbance stage and separately for the hurricane stage. 
In addition cross-sections were also plotted, in which the azimuth was the 
abscissa and standard pressure level was the ordinate. This was done for 
the distance ranging from 400 km to 800 km; this distance is well beyond 
most central overcast areas which generally do not exceed 300 km in radius, 

The result indicated that temperature distributions about the storms 
varied by about 2 - 30 K, which agrees with the findings of other investi- 
gators. There was a tendency for the western (forward) side of the 
hurricanes to be about 2°K warmer than the eastern side. Moreover, in 
Hurricane Camille, the hurricane stage was about 2°K warmer than the 
"disturbance'' stage in the western part of the storm. 

If these results are typical, and ^ the changes typically occur 
in sufficiently thick layers (at least 200 mb thick) , so that it might 
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be possible to detect the changes from satellite radiance measurements, 
and if the retrieval process can detect changes of 2-3°K in the presence 
of inevitable "noise" sources, the satellite data might be able to serve 
as a forecasting aid. To pursue this further, it will be necessary 
eventually to examine real satellite data over tropical disturbances 
and hurricanes together with airplane and other types of data. 
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Fig. 1. Vertical temperature profiles. Omaha temperature 
profile was used as a "first guess". Springfield temperature 
profile was used as the "true" temperature, in order to com- 
pute the "measured" radiance , and to judge the retrieval error, 
at Springfield. Ihe minimum-information retrieval profile was 
computed by using Omaha as the "first guess" and the simulated 
"measured" radiance at Springfield. 
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Fig. 2. Temperature errors, at Springfield, in the 
minimum-information retrievals as a function of the 
noise parameter The temperature errors were 

averaged over the pressure Intervals indicated on 
the curves. The first guess errors , for each layer, 
are indicated on the abscissa. 


S 
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g. 3|(a) . Coefficients F(k,j) used in the minimum-inf or- 
mation retrievals with Omaha temperatures as first guess 
Only 3 of the six spectral intervals, j, are shown. 

Curve (1) is for channel (1), v = 668.5 cm"^; curve 3, 

V == 695.3 cm“^; curve 6, ■o = 747.6 cm“^. (v = spectral 
frequency) . 






F(k,i) or G(k,j) 


Fig. 3(b). Coefficients for channel 6, v = 747.6 cm”^. Curve (3) 
is for F(k,6) used in the minimum-inf omiation retrievals; this 
is the same as curve (6) in Fig. 3(a). Curve 4 is for G(k,6), 
used in adjustment method (2); curve (A) is for Method (3); 
and curve (3) is for Method (4) . NOTE that all G(k,6) are . 
substantially smaller than F(k,6) , for k = 100, p = 1000 mbs. 
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Fig. 3(d) Coefficients for channel Ij v - 668.5 cm 

Curve (C) is for P(k,l) used in the minimim- information 
retrievals; this is the same as curve (1) in Fig. 3(a). 
Curve (A) is for G(k,l) for adjustment method (3); 
and curve (B) for Method (4). 





AVERASE ABSOLUTE TEMPERATURE ERROR *X 


g. 4. The average absolute temperature errors for 
indicated layers. On the ordinate, "900" means the 
layer from 1000 to 900 mbs; "800" means 1000 to 800 mbs; 

"700" means 1000 to 700 mbs; etc. The curves labelled 

(2) , (3), (4) refer to the "adjustment" retrieval methods (2) 

(3) , (4) . (see text) . 
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Fig. 5. Average absolute height errors at indicated 
pressure levels. The curves labelled (2) , (3) , 

(4) refer to the '’adjustment" retrieval methods (2), 
(3) , (4) . (see text) . 
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Fig. 6(a) and 6(b) . Temperature errors, averaged over indicated 
pressure layers, at '’stations” between Omaha and Springfield. 
Profile No. 2 means station at "Omaha +90 km”; No. 3 means 
"Omaha +180" km; etc.; No. 6 means "Omaha +450", or 90 km 
west of Springfield. Minimum-information method compared 
with adjustment method (4) . 










ABSTRACT 


This paper presents a method for using satellite measurements to 
interpolate vertical temperature soundings between radiosonde stations. 

The calculations presented show that especially in the 1000-800 mbs 
layer, where linear methods of temperature retrieval usually contain 
large errors, the proposed method reduces the errors substantially. 

The method finds a set of coefficients, which when multiplied by 
corresponding measured radiance quantities, yield zero temperature error 
at a radiosonde station. This derived set of coefficients is then ap- 
plied to satellite radiance measurements at places between the radiosonde 
stations. The computations show, for example, that the average absolute 
error in the layer 1000-800 mbs is only O.SKwhen the corresponding 
"minimum-information" method error was 2.9K. The method may be most 
applicable to measurements from geostationary satellites, but should also 
be applicable to measurements from polar orbiting satellites under certain 


conditions . 
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1. Introduction 


Indirect soundings of atmospheric temperature from satellite radiance 
observations are now obtained operationally from polar orbiting satellites 
(McMillin, et al. , 1973). In addition plans are going forward to mount a 
’’sounding" instrument on a geostationary satellite (Shenk, 1976). 

Two temperature retrieval methods have been xjsed in operational prac- 
tice. One is the so-called "minimum-information" method. The second is 
a regression method. Both are linear methods (Fritz et al. , 1972, 

Fleming and Smith, 1972) and seek to find the "best" set of coefficients 
F(k, j) in the following relationship: 


ABj.(k) - Z F(k, j)-AR(j) 

(1) 

2 


B^(k) - Bg(k) + AB^(k) 

(2) 

AR(j) “ R (j) - R (j) 
m s 

(3) 


In these equations 

B(k) “Planck function at the pressure level, k, at a reference fre- 
quency, usually taken at 700 cm when the 15 pm CO^ band is used in the 
satellite measurements . 

R(j) is a radiance at the frequency, j . 

The subscripts denote the following; 

”g" “ a guess (or initial estimate) 
m = the measured value 
r = the retrieved value. 

Because B(k) depends only on the tenqjerature , T(k) , it is easy to 
compute T(k) from the inverse of the Planck function, once B(k) is known. 

- 
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Unfortunately, the derived or retrieved temperatures often contain 
fairly large errors; in actual practice the RMS errors when compared with 
radiosondes are usually about 3K in the layer from 1000 to 800 mbs and 
also near the 200 mb level (Fritz, 1974; Weinreb and Fleming, 1974; 
vjerbowetski, 1975). Even with idealized, simulated "data", the average 
absolute error is usually about 3fC (Fleming and Smith, 1972). In part 
these errors in temperature retrievals are related to the shape of the 
vertical temperature profile (Fritz, 1969); the retrievals are not able 
to reproduce discontinuities well, and such discontinuities occur at the 
ground and at the tropopause. 

In this paper, a method is proposed which can improve the ten?>erature 
retrieval accuracy especially in the lower atmosphere and probably near 
the tropopause also. The results will be compared with computations 
derived from the minimum-information method. 

2. The minimum-information method 

In all temperature retrieval methods , satellite radiances are measured 
at a set of frequencies. We shall use the six CO 2 15 pm frequencies ob- 
served by the VTPR (Vertical Temperature Profile Radiometer) on NOAA 
satellites (McMillin et al. , 1973). 


TABLE 1 

Frequencies of channels used (cm ^) 


Channel No . 

1 

2 

3 

4 

5 

6 

Frequency 

668.5 

677.9 

695.3 

708.6 

725.5 

747.6 
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In the mlnimum-in formation method the values of F(k, j) in Eq. 1 are 
given by (Fritz, et al., 1972, Fleming. and Smith, 1972): 


5’(k. j) = 


(4) 


where -1 denotes the inverse matrix 

A = the matrix of 3) 

dx 

T 

A =■ the transpose of A 
2 ^ 

* the co-varxance of temperature between levels k in an ensemble 
of temperature profiles, but assumed constant. 




the varxance of noise in the ob4v:rvations ; instrument noise assumed. 


T(k, j) » the atmospheric transmittance from the level, k, to the top 
of the atmosphere, at the spectral frequency, j. 

X =» a height parameter, taken as proportional to 
In the minimum-information method, is taken as constant although 
in fact it varies with height. For the instrumental noise is usually 

assumed, although other factors, such as clouds, introduce "noise" into 
the observations. In this study = 100 was assumed, and = 0.50 for 
channel 1 and = 0.25 for the five remaining channels. [The transmit- 
tances, and temperature corrections for the transmittance were kindly 
supplied in a computer program by Dr. M. P. Weinreb (NESS/NOAA) ] . The 
program contains 100 levels for k, and six spectral intervals for j. 

The inverse matrix in Eq. (4) is so complex that, even for well behaved 
functions it is not possible to predict the values of F(k/ j) in advance. 


In addition, the values of F(k, j) are highly sensitive to values of 


N 


and selected. 


It is likely therefore that near any particular place and time, the 
set of F(k, j) may not be the "best". Furthermore, because of the 
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assuniptions adopted in the minimxim information method, the set of F(k, j) 
may not even be the "best” set for an ensemble. 

Even if AR(j) = 0 in Eq. (1), the retrieval B^(k) will in general not 
agree with the true value, B^(k) , at a level k, although B^(k) may oscil- 
late about the true values in some- unpredictable manner. In particular, 
as stated earlier, large deviations from the true values appear in the 
lower atmosphere and near the tropopause. 

It therefore, seemed reasonable to modify the set of E(k, j) in order 
to force B^(k) to equal the true values, B^(k) , for at least one sounding. 
After the new set of coefficients, G(k, j) had been selected, giving 
errorless retrievals for a known temperature profile at a radiosonde sta- 
tion, the set of G(k, j) could be used to interpolate between radiosonde 
stations, or extrapolate in time between the times of radiosonde observa- 
tions. This interpolation and extrapolation capability might-be especially 
useful when continual radiance observations become available from geosta- 
tionary satellites. The interpolation capability, between radiosonde 
sounding stations, might also be useful for radiances observed from polar 
orbiting satellites, since the field of view of the satellite radiometers 
i.s often much smaller than the distance between radiosonde stations. 
However, such interpolation and extrapolation capabilities would be most 

i 

useful in rapidly changing, small scale situations, provided the tempera- 
ture changes were large enough ; for this reason the application to the 
continual observations from geostationary satellites might be more pro- 
mising. Unfortunately such situations are often accompanied by complicated 
cloud fields, which may introduce errors into the "clear-colimm" radiances 
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required in most retrieval methods. Nevertheless, in this paper, it will 
be assumed that accurate "clear-column” radiances are available. 

3. Adjusting the coefficients; the adjustment method 
Let us find a set of coefficients, such that 


AB'^(k) = ABj.:(k) =» S G(k, j) *AR(j) 

j 

B’^(k) = B^.(k) » B^(k) + AB^Ck) ■ 


(5) 

( 6 ) 


The subscript, t, denotes the true value of the quantity. 

If we divide Eq. (6) by Eq. (2), 

AB^(k)/AB^(k) = [B fk) - B^(k)]/[B^(k) - B^(k)] = C(k) (7) 

^ t w g X 9 


The quantity C(k) will be a constant for a particular level, after the 
minimum-information retrieval determines the value of B^(k) . Substituting 
for AB^(k) from Eq. (1) , we obtain 


ABj.(k) = C(k)-AB^(k) = I C(k)-F(k, j)-AR(j) (8) 

j 


= 2 G(k, j) -ARCj) 


(9) 


G(k, j) = C(k)-F(k, j) 


( 10 ) 


4. Test of the Adjustment Methods 

To test the proposed method, several numerical testa were made. All 
the tests improved the minimum- inf ormation retrievals, especially in the 
lower atmosphere. One expects the satellite radiances to yield the most 
significant results where large temperature variations occur. Therefore, 
a tornado producing situation, with large horizontal temperature gradient 
was selected. Danielsen (1975, see p. 180) shows a vertical cross-section 



of potential temperature which contained a marked gradient between Omaha, 
Nebraska and Springfield, Illinois. These stations are about 550 km 
apart. From Danielsen's cross section, vertical temperature profiles were 
computed at five additional points located at every 90 km befcJ’een Omaha 
and Springfield. Thus seven soundings were available, including Omaha 
and Springfield. Above the top (about 160 mbs) of Danielsen’s diagram, 
an arbitrary climatological temperature sounding, typical of the lati- 
tude and time (April), was added up to the 0.01 mb level (about 80 km). 

(This climatological temperature add-on was also supplied by Dr. Weinreb) . 
For each of the seven locations the radiance was computed from the tempera- 
ture profiles from the radiative transfer aquation, 

R(j) - 2 B(k, j)*AT(k, j) + Bg(j)Xg(j) + NCj) (11) 

k 

where Bg(j) is the Planck function corresponding to the surface temperature 
and transmittance through the entire atmosphere at the fre- 

quency, j. The B(k, j) are computed from the temperatures, T(k) , from the 
Planck fmction formulai. The quantity N is a random nimber added to repre- 
sent instrumental noise. This was selected from a sub-routine which 
generated random numbers with a mean of zero and a standard deviation of 
0.5 for Channel 1, and 0.25 for the other channels. The experiment was 
also run with N = 0. The results were essentially the same. 

The radiances, R(j), in Eq. (11) were normalized to radiances at 700 cm 
so that in the minimum-information retrieval, quantities involving B(k) in 
Eq. (1) have all been scaled to 700 cm the temperatures can then be 
readily computed from the Planck function B(k,700) vhich are also the B's 
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With, these simulated ’’measured" radiances, » fi’o® Eq.(ll) , the 

minimum-information retrievals were obtained as follows; 

The Omaha somding was used as the ’’first-guess’’ or initial estimate. 

This supplied values of B (k) and of R (j) . With the aid of the Cciaha 

S S 

temperatures, and the t’s, a set of F(k,j) was computed from Eq.(4), With 
these values and the R^(j) for Springfield B^(k) was obtained from Eq3.(l) 
and (2) ; these are the retrieved values for Springfield. 

Using the Springfield temperatures as the truth, we compute C(k) from 
Eq.(7); and from G(k) and F(k,j) we compute 

G(k,j) from Eq. (10). 

These values of G(k,j) , when used with the Springfield radiances and Omaha 
temperatures as the first guess, yield a retrieval with zero error at all 
pressure levels at Springfield. 

Now these same values of G(k,j) and the same Omaha temperature sound- 
ing as first guess were used to obtain retrievals at the five locations 
between Omaha and Springfield, For comparison the retrievals were also 
obtained at those five locations with the minimum-information method, 

' I 

also using Omaha as the first guess. 

The results* are shown in Table 2. The column labelled "adjustment 
method" is the method involving G(k,j) . 


*The results shown here were not obtained by iteration; when the minimum- 
information solutions were obtained by iterations, the results were often 
poorer. 
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errors at all the "stations", for the 510 mb pressure level (one of the 
"k" levels was at 510 mb) . In all cases except the "Omaha +90 km" 
station, the height errors were smaller for the "adjustment" method. 
However, the height gradients, illustrated by iE(z) in Table 3, were 
fairly similar for the two methods- of retrieval. It would appear there- 
fore that if both satellite data and radiosonde data were used without 
adjustment to determine the height field, there would be a marked discon- 
tinuity near some radiosonde stations (e.g., Springfield) when compared 
with minimum-information retrievals. With the adjustment method however, 
the height field at nearly every station has been moved up to more nearly 
agree with the true heights. Of course this adjustment is forced to be 
exact at Springfield so that both sa tellite data and radiosonde data agree 
there. 

Finally, it should be pointed out, that a linear interpolation of the 
height field between Omaha and Springfield contains even larger height 
errors than those shoxra in Table 3; in this example, the variation of 
height between Omaha and Springfield was non-linear, and the satellite 
"data" yielded heights closer to the truth than linear interpolation of 
the height field did. 


TABLE 3 

Height Errors (m) for the 509 mb Pressure Surface 


Station 

Omaha 

+90 

km 

+180 

+270 

+360 

+450 

Springfield 

First Guess (Deviation 
from True) 

(Omaha First Guess) 

0 


-29 


-115 


-135 


-139 


-14] 


-142 

(Minimimi-Inf 0 rmation) 
AE(z) ; 

0 

+14 

14 

-49 

-35 

' ' 

-7 

-42 

+6 

-36 

+8 

-28 

0 

-28 

Adjustment ' 
AE(z) 

0 

+23 

23 

-39 

-16 

-5 

-21 

+9 

-12 

^'ll 

-1 

+1 

0 


AE(zj =■ the error in difference of heights between adjacent stations 
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5. Discussion and Conclusion 

We have demonstrated that if accurate radiances can be obtained from 
satellite measurements, then it is better to adjust the minimum-information 
coefficients, F(k,j), so that at least one satellite temperature retrieval 
sounding agrees exactly with its corresponding radiosonde temperature 
sounding, than to use the minimum-information retrievals themselves. This 
would apparently be especially valuable in the lower atmospheric layers 
where large temperature errors appear both in operational results and in 
other theoretical studies. For the two classes of cases we studied, the 
Gate ship soundings (not reported here) and the Omaha-Springfield tornado 
situation, the coefficients F(k,j) were too large by a factor of 3 or 
more near the 1000 mb level. This raises the question as to whether the 
assumptions, which enter the minimum- information method, always give too 
much weight to the coefficients in the lower layers, 

4 

However, the method proposed here, the Adjustment Method, has some 
problems too. For example, the factors C(k) in Eq.(7) may become un- 
stable if B (k)= B (k) . Substantial oscillations of G(k) were encoun- 
tered in this study too. But these were confined to narrow layers, and 
always oscillated in sign in these layers. This sign variation, of 
course, reduces the oscillation to small amplitude when temperature 
averages over layers are calculated. If it is necessary to reduce the 
temperature error at every level, k, then the G(k,j) can be averaged 
before the temperatures are calculated. This will no longer produce 
zero error at every level, k, at the adjustment station (Springfield 
in this Study), but the error will still be small, and the error over 
layers will still be nearly zero, because of the linearity of the opera- 
tions. Alternatively, additional methods may be developed to compute 
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G(k,j) such as limiting S G(k,j) to some constant value. Such methods 

j 

are being investigated further. 

The Adjustment Method discussed here suggests application to geo- 
stationary satellite radiances. For in that case frequent observations 
with relatively high spatial resolution will become available. Therefore 
observations would be available to study phenomena which are changing 
rapidly between the 12 hour radiosonde interval and also phenomena which 
are smaller in size than the distance between radiosonde stations. How- 
ever, clouds and water vapor will introduce errors in the "clear-column” 
radiances in these small time and space scale phenomena. The reduction 
of such radiance errors will present a challenging problem. 

Furthermore, it may also be possible to improve the temperature sound' 
ings from polar orbiting satellites with the proposed methods. The ref- 
erences already cited show that the largest errors occur near the earth's 
surface and also near the tropopause. In the tropics where the shape of 

the vertical temperature profile does not change much with time or space 

1 - ^ 

it might be desirable to use the adjustment method. For example, we 
might be able to take the Barbados sounding as a first guess, and re- 
trieve the Dakar temperature profile with the minimum-information method. 
Then adjust the coefficients F(k,j) so that the Dakar sounding is matched 
with no error. Finally, use- the new coefficients, G(k,j), to retrieve 
temperatures for Atlantic Ocean areas in between, and also into the 
Caribbean Sea area. Verification can be checked with the other Caribbean 
Sea area radiosonde stations and with the island stations near Dakar. 
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Care must of course be taken that the radiosonde pair selected are 
themselves on the same scale; for apparently radiosondes from different 
countries may not all be on the same scale. At first, it might be better 
to test with the Gate or Bomex ship data, to avoid discrepancies between 
the low level island temperatures -and temperatures over the sea. Such 
procedures may improve the operational satellite temperature retrievals 
both near the surface and near the tropopause. 
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APPENDIX B 

Methods of adjusting F(k,i) to G(k,j) (see p. 4) 
so that T^(k) - T^(k) = 0 at Springfield, Illinois 

Method (1) (see Appendix A) 

Change only F(k,j^^) 

Given: B^(k) = minimum-information retrieval at Springfield based on 

B (k) = first-guess Planck functions for 
S 

Omaha, Nebraska, the first-guess station, and 
AB (k) = SF(k,j).AR(j) 

^max ~ frequency for which [F(k,j)| has its maximum 


value at a given k. 

For the minimtim-information retrievals 

B^(k) = Bg(k) + AB^(k) (1) 

AB (k) = EF(k,j) -AR(j) (2) 

Let the true value of B be given by 

Bj.(k) = Bg + AB^(k) (3) 

subscript, t, means true value. 

Let AB (k) = EG(k,j) -ARCj) (4) 

^ j 

ABj.(k) - AB^(k) = [B^(k) - Bg(k)] - [B^(k) - B^Ck)] 

= B^Ck) - B^(k) (5) 

Therefore, from eqs . 2, 4, and 5 

ZG(k,j)-AR(j) - SF(k,j)*AR(j) = B. (k) - B^(k) 

A A LL 


I[G(k,j) - F(k.j)] .AR(j) = Bj.(k) - B^(k) 


or 


( 6 ) 
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Assume that G(k,j) = F(k,j) except for that value of j for which |F(k,j)I 
has a maximum value; this maximum value is defined separately for each 
value of k. Therefore, from eq. 6, 

= B,(k) - S^(k) 


or 


B. (k) - B (k) 

G(k,i ) = , r + F(k,j ) 

AE.(j ) ■'max 


'max 


max 


= C(k) + F(k,j ) 
^ ^ ’-^max 


Method (2) 

Change all F(k,.j) 

As in Method (1) 

AB^(k) = B^Ck) -B^(k) = ?F(k,j)*AR(j) 
r r s j 

AB fk) = B. (k) - B (k) = ZG(k,j).AR(j) 
t I- O J 

Dividing eq. (2) by eq. (1) 



Bj.(k) - Bg(k) 
Br(k) - Bg(k) 


]AB^(k) 


- C(k) • B^(k) 

once B (k) is obtained from the minimum-information retrieval, 
r 

From eqs . 1-4 

ZG(k,j)vAR(j) = C(k)IF(k,j)*AR(j) 

j . j 

= SC(k)-F(k,j)vAR(j) 
j 


( 1 ) 

( 2 ) 


(3) 

(4) 
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Therefore 

G(k,j) = C(k)-F(k.j) 

is a relationship which will force = T^. at Springfield, Illinois if 
the temperature profile at Omaha, Nebraska is the first-guess. 

Method 3 

2G(k,i) = 1 
and 

Change the two values of F(k,j) for which the absolute values of F(k,i) 
are the largest. 

As in Methods (1) and (2) 

= i:F(k,j) *AR(j) (1) 

" 3 

AB = EG(k,j).AR(j) (2) 

3 


Let 


ZG('k,j) = 1 


(3) 


and adjust only I 1 and I F(k, I ; is the frequency 

for which [F(k,j)[ is a maximum; frequency for which 

|F(k,j)[ has the second largest value at the level , k. 

For all other values of j 
G(k,j) ^ F(k,j) 

Therefore 

AB^(k) VF(k,j^^^).AR(j^^^) + F(k,j^^_^)-ARCj^^^_p 
4 

+ EF(k,j) *AR(j) 

3 


(4) 


(5) 


AB^(k) = + G(k,j^^^.^),AR(j^^^.^ 


+ EF(k,j).AR(j) 
3 ' 


( 6 ) 
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and from eq. (3) 

4 

+ G(k, ) + EF(k,j) - 1 


'max 


(7) 


E = sum over the four values for which j j 

^ 2 ^ : 


max 


max-1 


Eqs . (6) and (7) constitute two equations in the two unknowns, ^^^»^niax^ 

and G(k, j ) . 

-'max-1 

Hiltiply eq. (7) by AR(j ,) and subtract from eq. (6). 


Then, 


AB 




4 4 

+ SF(k,j)-AR(j) - [EF(k,j)]-AR(j^^_j_) 
J . J 






( 8 ) 


Once jjjjax^ computed from eq. (8) , we compute 




Method 4 

In Method 3 the denominator in eq. 8 may become small when 

AR(j ) - AR(j t). 

^■^max ■'max-1 

When this happens the coefficient G(k,j „) may become rather large, 

uIoX 

and the retrieval temperature errors may become large when Che coefficients 
are applied to independent satellite ''measurements”. 
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To remedy this defect a procedure similar to the one for Method 3 

can be used. But now select R(i ) and in addition one other R(j) [not 

max 

necessarily R(i t)]* Select the second R(j) such that 
maX“l 


AR(i ) - AR(j) is a maximinn. 


The procedure is as follows : 

Find F(k,i ); this defines j and therefore AR(j ) . Now 
max max max 

test each AR(j) to see which one makes 
I AR(j^^^) - AR(j) I a maximum. 

For six measured channels, the newly selected "j" could take on 

any one of five values. After the second "j" has been selected, 

designated j ’ , proceed as in Method 3 except that j ' replaces 

i , . Then 
“^max-l 


G(k,j^^) = {ABj.(k) - AR(j') - IF(k,j)-ARCj) 


+ [2F(k,j)]-AR(j’)}/{AR(j^^) - AR(j*)}; j ¥ 

i‘ j 


f 


